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The BABAR and Belle collaborations have recently found evidence for mixing within the D meson system. We 
present some of the mixing search techniques used by BABAR and their status as of the beginning of the summer 
2007. These have culminated in a measurement in the Kn decay final state of the D that is inconsistent with 
the no-mixing hypothesis with a significance of 3.9 standard deviations. 
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1. Introduction 

Mixing among the lightest neutral mesons of each 
flavor has traditionally provided important informa- 
tion on the electroweak interactions, the CKM matrix, 
and the possible virtual constituents that can lead to 
mixing. Among the long-lived mesons, the D meson 
system exhibits the smallest mixing phenomena. The 
B-factories have now accumulated sufficient luminos- 
ity to observe mixing in the D system and we can 
expect to see more detailed results as more luminos- 
ity is accumulated and additional channels sensitive to 
mixing are analyzed. The B-factories produce about 
1.3 million Charm events per fb^^ of integrated lumi- 
nosity accumulated. The BaBAR integrated luminos- 
ity of about 384 fb~^ used for the evidence for mixing 
result we will present corresponds to about 500 mil- 
lion charm events produced. The present BaBAR inte- 
grated luminosity is approximately 500 fb^^. BaBAR 
is a high acceptance general-purpose detector provid- 
ing excellent tracking, vertexing, particle ID, and neu- 
trals detection. All of these capabilities are crucial for 
making the difScult mixing measurement. 



2. Mixing Measureables for the D 
System 

The propagation eigenstates, including the elec- 
troweak interactions for the D mesons are given by: 

\D,,,)=p\D^)±q\D"), \pf + \qf = l. (1) 

Propagation parameters that determine the time- 
evolution for the two states are given by: 

r = i(ri+r2), aa/ = Mi-M2, Ar = (ri-r2); 

(2) 

with the observable oscillations determined by the 
scaled parameters 



AM 



_ Ar 



(3) 



In the case of CP conservation the two D eigen- 
states are the CP even and odd combinations. We 



will choose Di to be the CP even state. The sign 
choice for the mass and width difference varies among 
papers, we use the choice above. 

Assuming CP conservation, small mixing parame- 
ters, and an initial state tagged as a we can write 
the time dependence to first order in x and y: 



D{t) = I +D°{-y-ix)-t 



-(r/2+im)t 



(4) 



Projecting this onto a final state / gives to first order 
the amplitude for finding /: 



Af + Af{-y - ix)^t 



-{T/2+vm)t 



(5) 



This leads to a number of ways to measure the effect 
of mixing, for example: 

(1) Wrong sign semileptonic decays. Here Af is zero 
and we measure directly the quantity, after integrating 
over decay times: 



i?Af = (.t2+2/')/2. 



(6) 



Limits using this measurement, however, are not 
yet sensitive enough to get down to the 10~* level for 
Rm- Using 334 fb~^ of data, electron decays only, 
and a double tag technique, BABAR measures Rm = 
0.4x 10~^, with a 68% confidence interval (—5.6, 7.4) x 
10-4[i]. 

(2) Cabibbo favored, right sign (RS) hadronic de- 
cays (for example K~t:'^). These are used to measure 
the average lifetime, with the correction from the term 
involving x and y usually ignored (provides a correc- 
tion of 0(10-3)). 

(3) Singly suppressed decays (for example K^K~ 
or TT+TT^). In this case tagging the initial state isn't 
necessary. For CP even final states: Af — Af. This 
provides the most direct way to measure y. With tag- 
ging we can also check for CP violation, by looking at 
the value of y for each tag type. BABAR will be updat- 
ing this measurement with the full statistics later this 
year. The initial BABMi measurement was based on 
91 fb~^ and gave the result y = 0.8%, with statistical 
and systematic errors each about 0.4% [2], measure- 
ment [3]. 
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4) Doubly suppressed and mixed, wrong sign (WS) 
decays (for example K~^Tr~ . Mixing leads to an ex- 
ponential term multiplied by both a linear and a 
quadratic term in t. The quadratic term has a uni- 
versal form depending on Rm- For any point in the 
decay phase space the decay rate is given by 

(^\Aff + \Af\ \Af\ yTt + \Af\^Rjlf-^ e"". 

(7) 

Here y' ^ y cos S — x sin 6, where (5 is a strong phase 
difference between the Cabibbo favored and doubly 
suppressed amplitudes. For the K'^n~ decay there is 
just the one phase and the ratio of \Af\^ to is 
defined to be Rd- For multibody decays the strong 
phase varies over the phase space and the term pro- 
portional to t will involve a sum with different phases 
if we add all events in a given channel. 

BaBAR has analyzed the decay channel K~^'K~n°, 
with a mass cut that selects mostly K^p~ decays, 
the largest channel for the Cabibbo allowed ampli- 
tude arising from mixing. Based on 230 fb~^, BaBMl 
measures [4] 

ay' = (-1.2tO;|±0.2)% 
Rm = {Omstoil ± 0.004)%. (8) 

The parameter a allows for the phase variation over 
the region summed over. A fit to the full Dalitz plot 
would allow more events to be used in the mixing 
study. This, however, requires a model for all the 
resonant and smooth components that contribute to 
the given channel, which may introduce uncertainties. 
BaBAR is working on such a fit, which will be based 
on approximately 1500 signal events. 

Another important 3-body channel is the Ksn~^n~ 
decay channel. Analysis of this channel was pioneered 
by CLEO [5]. It contains: CP-even, CP-odd, and 
mixed-CP resonances. Now one must correctly model 
the relative amounts of CP-odd and CP-even contri- 
butions (including smooth components) to get the cor- 
rect lifetime difference. This channel also provides the 
possibility to directly measure x. BABAR is working on 
this channel; Belle has published their results [6]. 

In the Standard Model y and x are mainly due 
to long-distance effects. They may be comparable in 
value but this depends on physics that is difficult to 
model. Long-distance effects control how complete the 
SU(3) cancellation is, which would make both param- 
eters vanish in the symmetry limit. The exact values 
therefore depends on SU(3) violations in matrix ele- 
ments and phase space. Also, the sign oi x/y provides 
an important measurement. One might expect the x 
and y parameters to be in the range 0(10^"' to 10^^). 
Thus the present data are consistent with the Stan- 
dard Model. Searches for CP violation are important 
goals of the B-factories, since observation at a non- 
negligible level would signify new physics. 



We will turn now to the strongest evidence for D- 
mixing from BaBAR, using the Ktt final state. This 
result has recently been published [7]. 



3. Analysis of the Ktt channel 

We study the right-sign (RS), Cabibbo- favored 
(CF) decay D° K~tt~^ [8] and the wrong-sign 
(WS) decay D° K+n-. The latter can be pro- 
duced via the doubly Cabibbo-suppressed (DCS) de- 
cay D° — > K~^TT~ or via mixing followed by a CF decay 
_D0 ^ ^0 ^ X+TT-. The DCS decay has a small rate 
Rd of order tan'' Oq ~ 0.3% relative to the CF decay 
with 9c the Cabibbo angle. We tag the D'^ at pro- 
duction using the decay D*+ T^t^^ where the 7r+ 
is referred to as the "slow pion" . In RS decays the 7r+ 
and kaon have opposite charges, while in WS decays 
the charges arc the same. The time dependence of the 
WS decay rate is used to separate the contributions 
of DCS decays from D^-D'^ mixing. 

We study both CP-conserving and CP-violating 
cases. For the CP-conserving case, we fit for the pa- 
rameters Rd, x'"^ 1 and y' . To search for CP violation, 
we apply Eq. (7) to the and D° samples separately, 
fitting for the parameters {P^, x' , y'^} for 1)° (-|-) 
decays and {—) decays. 

We select Z?" candidates by pairing oppositely- 
charged tracks with a K^-k^ invariant mass rriK-K be- 
tween 1.81 and 1.92GeV/c^. We require the 7r+ to 
have a momcntiim in the laboratory frame greater 
than O.lGeV/c and in the e+e^ center-of-mass (CM) 
frame below 0.45 GeV/c. 

To obtain the proper decay time t and its error at 
for each Z?'^ candidate, we refit the and tt^ tracks, 
constraining them to originate from a common vertex. 
We also require the D° and 7r+ to originate from a 
common vertex, constrained by the position and size 
of the e+e" interaction region. The vertical RMS 
size of each beam is typically 6 /xm. We require the 
probability of the vertex-constrained combined fit 
P(x^) to be at least 0.1%, and the vriK-K-n^ —rn,K-K mass 
difference Am to satisfy 0.14 < Am < 0.16GeV/c^. 

To remove candidates from P-meson decays and 
to reduce combinatorial backgrounds, we require each 
to have a momentum in the CM frame greater 
than 2.5 GeV/c. We require —2 < t < 4ps and 
Gt < 0.5 ps (the most probable value of at for signal 
events is 0.16 ps). For _D*+ candidates sharing one or 
more tracks with other _D*+ candidates, we retain only 
the candidate with the highest P(x^). After applying 
all criteria, we keep approximately 1,229,000 RS and 
64,000 WS D° and D" candidates. 

The mixing parameters are determined in an un- 
binned, extended maximum-likelihood fit to the RS 
and WS data samples over the four obscrvables mK-m 
Am, t, and at- The fit is performed in several stages. 
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First, RS and WS signal and background shape pa- 
rameters are determined from a fit to ttikw and Am, 
and are not varied in subsequent fits. Next, the D° 
proper-time resolution function and lifetime are deter- 
mined in a fit to the RS data using ryiK-K and Am to 
separate the signal and background components. We 
fit to the WS data sample using three different mod- 
els. The first model assumes both CP conservation 
and the absence of mixing. The second model allows 
for mixing, but assumes no CP violation. The third 
model allows for both mixing and CP violation. 

The RS and WS {mxTr, Am} distributions arc de- 
scribed by four components: signal, random 7r+, mis- 
reconstructed and combinatorial background. The 
signal component has a characteristic peak in both 
rriKTr and Am. The random 7r+ component models re- 
constructed D° decays combined with a random slow 
pion and has the same shape in niKn as signal events, 
but does not peak in Am. Misreconstructed 1?° events 
have one or more of the decay products either not 
reconstructed or reconstructed with the wrong parti- 
cle hypothesis. They peak in Am, but not in mx-n- 
For RS events, most of these are semileptonic de- 
cays. For WS events, the main contribution is RS 
D'^ K~Tr~^ decays where the K~ and the 7r+ are 
misidentified as tt^ and if respectively. Combinato- 
rial background events are those not described by the 
above components; they do not exhibit any peaking 
structure in m^Tr or Am. 

The functional forms of the probability density 
functions (PDFs) for the signal and background com- 
ponents are chosen based on studies of Monte Carlo 
(MC) samples. However, all parameters are deter- 
mined from two-dimensional likelihood fits to data 
over the full niKn and Am region. 

We fit the RS and WS data samples simultane- 
ously with shape parameters describing the signal and 
random 7r+ components shared between the two data 
samples. We find 1, 141, 500 ± 1, 200 RS signal events 
and 4, 030 ± 90 WS signal events. The dominant 
background component is the random Tr^ background. 
Projections of the WS data and fit are shown in Fig. 1. 




• Data 

□ WS Signal 
] Random ji, 
iMisrecon. D° - 

3 Combinatorial 



0.15 0.16 

Am (GeV/c^ 



Figure 1: a) rriKTv for wrong-sign (WS) candidates with 
0.1445 < Am < 0.1465 GoV/c^ and b) Am for WS candi- 
dates with 1.843 < mKn < 1-883 GeV/c^ The fitted PDFs 
are overlaid. 



The measured proper-time distribution for the RS 
signal is described by an exponential function con- 
volved with a resolution function whose parameters 
are determined by the fit to the data. The resolution 
function is the sum of three Gaussians with widths 
proportional to the estimated event-by-event proper- 
time imcertainty at. The random tt^ background is 
described by the same proper-time distribution as sig- 
nal events, since the slow pion has little weight in the 
vertex fit. The proper-time distribution of the com- 
binatorial background is described by a sum of two 
Gaussians. one of which has a power-law tail to ac- 
count for a small long-lived component. The combi- 
natorial background and real decays have differ- 
ent (T/ distributions, as determined from data using a 
background-subtraction technique based on the fit to 
niKn and Am. 

The fit to the RS proper-time distribution is per- 
formed over all events in the full tukv and Am re- 
gion. The PDFs for signal and background in rriKx 
and Am are used in the proper-time fit with all pa- 
rameters fixed to their previously determined values. 
The fitted D° lifetime is found to be consistent with 
the world-average lifetime [9] . 

The measured proper-time distribution for the WS 
signal is modeled by Eq. (7) convolved with the res- 
olution function determined in the RS proper-time 
fit. The random 7r+ and misreconstructed back- 
grounds are described by the RS signal proper-time 
distribution since they are real D° decays. The 
proper-time distribution for WS data is shown in 
Fig. 2. The fit results with and without mixing are 
shown as the overlaid curves. 

The fit with mixing provides a substantially bet- 
ter description of the data than the fit with no mix- 
ing. The significance of the mixing signal is evaluated 
based on the change in negative log likelihood with 
respect to the minimum. Figure 3 shows confidence- 
level (CL) contours calculated from the change in log 
likelihood (— 2Aln/Z) in two dimensions {x''^ and y') 
with systematic uncertainties included. The likeli- 
hood maximum is at the unphysical value of x''^ = 
-2.2 X 10-4 and y' = 9.7 x IQ-^. The value of 
— 2Aln£ at the most likely point in the physically al- 
lowed region {x'^ = and y' = 6.4 x 10"^) is 0.7 units. 
The value of — 2Aln>C for no-mixing is 23.9 miits. In- 
cluding the systematic uncertainties, this corresponds 
to a significance equivalent to 3.9 standard deviations 
(1 — CL = 1 X lO-"*^) and thus constitutes evidence for 
mixing. The fitted values of the mixing parameters 
and i?D are listed in Table I. The correlation coefla- 
cient between the x' and y' parameters is —0.94. 

Allowing for the possibility of CP violation, we 

calculate the values of Rn = R^R^ and Au = 
{R^ - R^)/{R^ + R^) listed in Table I, from the fit- 
ted values. The best fit points (x'^^, y'^) shown in 
Table I are more than three standard deviations away 
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Table I Results from the different fits. The first uncer- 
tainty listed is statistical and the second systematic. 



Figure 2: a) Projections of the proper-time distribution of 
combined and WS candidates and fit result inte- 
grated over the signal region 1.843 < rriKn < 1.883 GeV/c^ 
and 0.1445 < Am < 0.1465 GeV/c^ The resuh of the fit 
allowing (not allowing) mixing but not CP violation is 
overlaid as a solid (dashed) curve, b) The points repre- 
sent the difference between the data and the no-mixing 
fit. The solid curve shows the difference between fits with 
and without mixing. 




-0.5 0.0 0.5 

x'^/10-^ 



Figure 3: The central value (point) and confidence-level 
(CL) contours for 1 - CL = 0.317 (la), 4.55 x 10"^ (2a), 
2.70 X 10'=^ (3a), 6.33 x 10"'^ (4a) and 5.73 x 10"'^ (5a), 
calculated from the change in the value of — 21n£ com- 
pared with its value at the minimum. Systematic uncer- 
tainties are included. The no-mixing point is shown as a 
plus sign {+). 



from the no- mixing hypothesis. The shapes of the 

2_j- I 

(x' , y'*) CL contours are similar to those shown in 
Fig. 3. All cross checks indicate that the close agree- 
ment between the separate and D° fit results is 
coincidental. 

As a cross-check of the mixing signal, we perform 

independent {rriK-K, Am} fits with no shared param- 
eters for intervals in proper time selected to have ap- 



Fit type 


Parameter Fit Results (/lO"^) 


No CP viol. 


or mixing J?d 


3.53 ±0.08 ± 0.04 


No CP 
violation 


y' 


3.03 ±0.16 ±0.10 
-0.22 ±0.30 ±0.21 
9.7 ± 4.4 ± 3.1 




Ro 


3.03 ±0.16 ±0.10 


CP 


Ao 


-21 ± 52 ± 15 


violation 


x"+ 


-0.24 ±0.43 ±0.30 


allowed 


y'+ 


9.8 ± 6.4 ± 4.5 






-0.20 ±0.41 ± 0.29 




y'- 


9.6 ± 6.1 ± 4.3 



proximately equal numbers of RS candidates. The 
fitted WS branching fractions are shown in Fig. 4 and 
are seen to increase with time. The slope is consistent 
with the measured mixing parameters and inconsis- 
tent with the no-mixing hypothesis. 



0.45 



0.4 



5 
EC 



0.35 



0.3 

-2 ' ' ' ' -1 ' ' ' 1 2 

t(ps) 

Figure 4: The WS branching fractions from indepen- 
dent {rriKTv, Am} fits to slices in measured proper time 
(points). The dashed line shows the expected wrong-sign 
rate as determined from the mixing fit shown in Fig. 2. 
The with respect to expectation from the mixing fit is 
1.5; for the no-mixing hypothesis (a constant WS rate), 
the is 24.0. 

We validated the fitting procedure on simulated 
data samples using both MC samples with the full 
detector simulation and large parametrized MC sam- 
ples. In all cases we found the fit to be unbiased. As 
a further cross-check, we performed a fit to the RS 
data proper-time distribution allowing for mixing in 
the signal component; the fitted values of the mixing 
parameters are consistent with no mixing. 

In evaluating systematic uncertainties in i?D and 
the mixing parameters we considered variations in the 
fit model and in the selection criteria. We also consid- 
ered alternative forms of the ttikt!, Am, proper time, 
and Gt PDFs. We varied the t and at requirements. In 
addition, we considered variations that keep or reject 
all D*+ candidates sharing tracks with other candi- 
dates. 
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For each source of systematic error, we compute the 
significance sf = 2t[lnC{x'^, y') - ln/:(.T'^, ?/;)]/2.3, 
where (x'^, y') arc tlic parameters obtained from the 
standard fit, {x'^,y'-) the parameters from the fit in- 
cluding the i"' systematic variation, and C the likc- 
hhood of the standard fit. The factor 2.3 is the 68% 
confidence level for 2 degrees of freedom. To estimate 
the significance of our results in [x'"^ , y'), wc reduce 
— 2Aln£ by a factor of 1 + Ssf — 1.3 to account for 
systematic errors. The largest contribution to this fac- 
tor, 0.06, is due to uncertainty in modeling the long 
decay time component from other D decays in the sig- 
nal region. The second largest component, 0.05, is due 
to the presence of a non-zero mean in the proper time 
signal resolution PDF. The mean value is determined 
in the RS proper time fit to be 3.6 fs and is due to small 
misalignments in the detector. The error of 15 x 10~^ 
on Ad is primarily due to uncertainties in modeling 
the differences between K'^ and K~ absorption in the 
detector. 

In conclusion we summarize the BABAR evidence 
for D'^-D'^ mixing. Our result is inconsistent with 
the no-mixing hypothesis at a significance of 3.9 stan- 
dard deviations. We measure y' = [9.7 ±4.4 (stat.) ± 
3.1 (syst.)] X 10~^, while x'"^ is consistent with zero. 



We find no evidence for CP violation and measure Ru 
to be [0.303 ± 0.016 (stat.) ± 0.010 (syst.)]%. The re- 
sult is consistent with Standard Model estimates for 
mixing. 
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